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Abstract 

With the global color symmetry model (GCM) at finite chemical potential, the 
density dependence of the bag constant, the total energy and the radius of a nucleon 
in nuclear matter is investigated. A relation between the nuclear matter density and 
the chemical potential with the action of QCD being taken into account is obtained. 
A maximal nuclear matter density for the existence of the bag with three quarks 
confined within is given. The calculated results indicate that, before the maximal 
density is reached, the bag constant and the total energy of a nucleon decrease, 
and the radius of a nucleon increases slowly, with the increasing of the nuclear 

*Work supported by the National Natural Science Foundation of China, the Major State Basic Re- 
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matter density. As the maximal nuclear matter density is reached, the mass of the 
nucleon vanishes and the radius becomes infinite suddenly. It manifests that a phase 
transition from nucleons to quarks takes place. 

PACS Numbers: 24.85.+p, ll.lO.Wx, 12.39.Ba, 14.20.Dh 

Keywords: Effective field theory of QCD, Bag model, Nucleon, Density dependence 
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1 Introduction 



It is well known that nucleons bound in nuclear medium alter their properties from those 
in free space. For instance, an effective mass (< Mq) should be introduced to the nucleon 
to describe the property of nuclear matter [T], and a swell hypothesis should be made to de- 
scribe the EMC effect well (see for example Ref. 0). Although how much those properties 
change is of fundamental interest in nuclear physics, it is still not clear now. Then, in the 
relativistic mean field theory (RMF) based on the a — u model, the effective nucleon mass 
in nuclear medium is determined with iteration on the fields of nucleon and mesons jTj. 
However, this model is valid only at the hadronic level, at which a nucleon is described as 
a point-like particle. Because of the importance of the substructure of nucleon observed in 
deep inelastic scattering experiments and predicted by quantum chromodynamics (QCD), 
it is imperative to employ the quark and gluon degrees of freedom to describe nuclear phe- 
nomena. Since it is now very difficult to make quantitative predictions with QCD at low 
and intermediate energy region, various phenomenological models based on the QCD as- 
sumptions, such as the bag modelsjSl Ej, quark-meson coupling (QMC) model[5 , and so 
on, have been developed. In bag models and the QMC model, the size and mass of a 
nucleon are represented by the radius and energy of the bag with quarks being confined, 
respectively. In order to take the modification of the volume energy on the mass of a nu- 
cleon into account, the bag constant is exploited in the models. To take into account the 
nuclear matter medium effect in the QMC model, a phenomenological density dependence 
of the bag constant in the medium has been introduced HI El 13 El- In the global color 
symmetry model (GCM)[TTl[T21[inilIll of QCD, even though the bag constant, the radius 
and the mass of a nucleon in nuclear matter have been evaluated ^H] in a consistent way, 
a Fermi gas approximation was taken for the relation between the nuclear matter density 
and the chemical potential. The interaction among the nucleons which is believed to be 
the residual of the quark-quark interactions has not yet been taken into account. Then it 
is still necessary to investigate the nuclear matter density dependence of the radius, the 
mass and the bag constant of a nucleon sophisticatedly. 

It has been shown that the global color symmetry model (GCM) [TT| IT^ IT^ ITij is 
a quite successful effective field theory of QCD in describing hadron properties in free 
space(i.e., at temperature T = 0, chemical potential fi = 0). With the global color 
symmetry model at finite chemical potential /i, we will determine the relation between 
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the nuclear matter density and the chemical potential in the GCM and reevaluate the 
nuclear matter density dependence of the bag constant, the total energy and the radius 
of a nucleon in this paper. 

The paper is organized as follows. In Section 2 we describe the formalism of the GCM 
model at finite chemical potential /i and the relation between the chemical potential and 
the baryon density in nuclear matter. In Section 3 we represent the calculation and the 
obtained results of the bag constant, the bag radius and the bag energy as functions of 
the nuclear matter density. In Section 4, a brief summary and some remarks are given. 



2 Formalism 



The starting point of the global color symmetry model (GCM) is the action in Euclidean 
matricPT] 
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S = j d'^xd'^y 



q{x) (7 ■ 9 + m) 5(x - y)q{y) - ^j^{x)D^,^{x - y)il{y) 



where = q{x)^'~)^q{x) is the quark current, D^^{x — y) is an effective two-point gluon 
propagator, m is the current quark mass, g is the quark-gluon coupling constant. Taking 
the effective gluon propagator to be diagonal, i.e., D^^{x — y) = 6fj,uD{x — y) and applying 
the transformation of Fierz reordering to the quark fields, one can rewrite the current- 
current term as 
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2 J dVyj;{x)D{x-y)fM = J dVyf{x,y)D{x-y)f{y,x), 

where j^{x,y) = q{x)A.^q{y) with being the direct products of Lorentz, flavor and 
color matrices which produce the scalar, vector, pseudoscalar and axial vector terms. It 
is obvious that such a j^{x,y) is a bilical current. With two flavors u and d of quarks 
being taken into account, each A is either isoscalar or isovector. The color matrices in- 
volved in the Fierz transformation contains color-singlet and color-octet terms. Taking 
the bosonization procedure one can transfer the bilocal quark current structure into aux- 
iliary Bose-fields carrying the quantum number 9. The action of the GCM in free space 
(i.e., at the chemical potential = 0) for the zero-mass quark can then be rewritten |1 1] 
in the Euclidean space as 

S{B) = J d'xd^qix) [7 ■ d5{x -y) + A'B\x, y)] q{y) + j d'xd'y?^^^^^^ , 



where B^{x,y) is the bilocal Bose-field. As a consequence, the generating functional is 
given as 

Z[rJ,r]]=^^ I DgDgD5''e[-^(^'')+/'''"(^''+^^)l, 

where rj and rj are the quark sources. To extend the GCM to at finite nuclear matter 
density (with finite chemical potential fi) , one should, in view of the statistical mechanics, 
take the partition function of the canonical (quark) ensemble into that of the grand 
ensemble with quarks and hadrons that are the solitons collecting quarks jT2]. The quark 
field should be transformed under a constraint on the baryon number through the chemical 
potential /i 

q{x) > q'{x) = ef'^'^qix) . (1) 

After some derivation, we have the action of the GCM in nuclear matter 

5(5^/i) = Jd^xd^yq'{x)[{j-d-fx-fi)5{x-y) + e'''''A^B%x,y)e-^'y'] q'{y) 



+ J d^xd^y 



B%x,y)B%y,x) 



2gWix-y) 
and the generating functional is given as 

r],7]]= J DqDqDB^e^-^^^'^ . (2) 

After integrating the quark fields, we obtain 

5(5^/i) = -Trln[(7-9-/i74)(5(x-i/) + e'^"^A^5V^f*] + / rfWy^^^^J^^^^ . (3) 

J 2g D[x y) 

Generally, the bilocal field B^{x,y) can be written asfH] 

B\x, y) = Bl{x, y) + Y: r^(x, i/)0^(^) , (4) 

where BQ^x^y) =Bq{x — y) = B^{x — y) is the vacuum configuration of the bilocal field. 
SjFf (x, (^^) correspond to the fields which can be interpreted as effective meson 
fields. In the lowest order approximation with only the Goldstone bosons being taken into 
account, the 0f includes tt and a mesons. The corresponding width of the fluctuations are 
approximately the same as the vacuum conflguration|14j . i.e., = Bq. Since the bilocal 
fleld arises from the bilocal current of quarks, the internal q-q structure of the mesons 
can be well described in the GCM. The vacuum conflguration can be determined by the 
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saddle-point condition = 0. Then an equation of the translational invariant quark 
self-energy T,{q,fi) is obtained as a truncated Dyson- Schwinger equation 

J [271)^ 2 I'j ■ q- /i74 + E(g, /i) 2 

With q^ = (g, g4 + i/i) being introduced, the above equation can be rewritten as 

f d'^q 1°" 1 

= / 77r^9^D{p-q)--f^- ^——7^-. (5) 

J (27r)* 2 z7-g + L(g) 2 

Considering the fact that the inclusion of the chemical potential breaks the 0(4) symmetry 
in the four-dimensioanl space, one should rewrite the decomposition of the self-energy E 
as 

S(p, /i) = i[A{p, /i) - 1]7 ■ p + i[C{p, /i) - 1]74(P4 + «/i) + B{p, /i) , (6) 
where B{p,fi) is the counterpart of the vacuum configuration in the momentum space. 



I.e., 



i2ny 

It requires then that the vacuum configuration of a and tt should satisfy the restriction 
cr^ + vf^ = 1. Combining Eqs. (5) and (6), one can obtain the equations to determine the 
A{p), B{p) and C{p) as follows 

[A{p) - l]p = - / ——g D[p - q) ^^^^^ „ ^ ^o/;,x;,2 , m/^N ' (7a) 



3 7 (27r)4" ^' A^{q)q' + C\q)qi + B\q) 

\n(~\ 11 ~2 8 /■ 2n^ A C{q)q4P4 

= y y (2^0^^ - ^^IWTcWT^^ • ^'^^ 

Basing on the solution of the Dyson- Schwinger equations, one can determine the 
bilocal field, and fix further the GCM action. 

With a nontopological-soliton ansatz, the action of the bag models at finite chemical 
potential in terms of the GCM can be given by extending the formalism proposed in 
Refs.m] and as 

SB = qj{l-d--i4ii-a[a{x)~m{x)-T-i^]}qj + S{a,'K,ij), (j = 1,2,3), (8) 
where S{a,7i,fi) includes a and vr mesons and reads 

S{a, vr, /i) =1 l^my + ^-f{d,af + ^m^vr^ + fl[d,^f-V{a, vr)]rfS + ■ ■ • (9) 



with 



V{a,Tx) 



-12 



In 



A^{p)f + C^{p)pl + B^{p) 
{a^ + ^^-l)B\p) \ 



A\p)p' + C\p)pI + B\p) 
and the quark meson couphng constant a is given as 

d'^p 



a[x) 



-Mp)e 



ip-x 



It is evident that such a quark meson couphng constant is just the vacuum configuration 
of the bilocal field and is independent of the meson fields. 

Analyzing the stationary property of the bag and differentiating Eq. (8), one has the 
equations of motion for the quarks and mesons 



{7 ■ a - 74yu - a[a{x) - i7f{x) ■ tj^]} Qj = , 



(10) 



dS, 



B 



da{x) 



(11) 



dS 



B 



dn{x) 



0. 



(12) 



The quark field and a, tt meson fields in symmetric nuclear matter can be determined by 
solving the Eqs. (10-12) self-consistently. As a consequence, the corresponding energies 
can be obtained. It is apparent that the meson fields corresponding to the vacuum con- 
figuration can be simply taken as a = 1, vr = due to the restriction vr^ -|- o"^ = 1. In light 
of the nontopological-soliton ansatzjUE], one can approximately take the meson fields 
inside a bag (i.e., in a nucleon) as a = and vr = 0. For the quarks in a bag, Eq. (10) 
can thus be rewritten approximately as 



[7-9-74/i]gj(x) = 0. 



(13) 



The lowest total energy of a single quark with respect to the radius R of the bag is given 
as 



^0 - 
R ' 



(14) 
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where uq = 2.04. For a free nucleon (i.e., with fi = 0), this result is consistent with those 
obtained in RefsjHl El EI- And the bag constant B is obtained as 



i3 = 1^(0,0,) -1^(1,0) = 12 J- 



(2^ 



In 



A\p)f+C^ip)pl+B\p) 
A^{p)f+C'^{p)pl 
B\p) \ 



A^{p)f+C\p)pl+B^{p) 



(15) 



With the correction from the motion of center-of-mass, the zero-point effect and the color- 
electronic and color-magnetic interactions being taken into account, the total energy of a 
bag is given as 

E = 3e, {R) + \.R^B - | = "^^ ~ ^° + \.R^B , (16) 

where Zq/R denotes the corrections of the motion of center-of-mass, zero-point energy 
and other effects. 

Just as the same as that in Ref.|ll|, the bag is identified as a nucleon in the present 
work. It satisfies the equilibrium condition 

dE{R) 



dR 







From this condition, we get 



where a = 3{u!o — /i) — Zq. As a consequence, Eq. (16) can be rewritten as 

4a /47ri3\^/^ 



E 
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It is apparent that, with the solutions of Dyson- Schwinger equations (Eqs. (7a), (7b) 
and (7c)) being taken as the input for Eqs. (15), (17) and (18), the properties of nucleons 
(i.e., bags) in nuclear matter can be obtained. 

In the practical calculation, since the knowledge about the exact behavior of and 
D{p — q) in low energy region is still lacking, one has to take some approximations or 
phenomenological form to solve the Dyson-Schwinger equations. For simplicity, we adopt 
the infrared dominative formjini EI 

g^D{p-q) = ^r]^6ip-q), (19) 
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where rj is an energy-scale parameter and can be fixed by experimental data of hadrons. 
Although this form does not include the contribution from the ultraviolet energy region, 
it maintains the main property of QCD in the low energy region. With Eqs. (7a), (7b), 
(7c) and (19), one has 

-,2 



A{p) = C{p) 



A{p) = C{p) = - 



1+14 



B{p) = {v'-4p'y/\ 

2\ 1/2' 

B{p) = 0, 



27] 



p2 



for i?e(p2) < ^, (20a) 

2 

for i?e(p2) > ^, (206) 



which describes the phase where chiral symmetry is spontaneously broken and the dressed- 
quarks are confined |12[ ITTl ITHl I19| I20j. Meanwhile one has also the Wigner solution 



A(P) = C(p) 



p 



2\ 1/2' 



B{p) = 0, 



(21) 



which characterizes a phase in which chiral symmetry is not broken and the dressed-quarks 
are not confined. 

In order to investigate the dependence of nucleon properties on the nuclear matter 
density p explicitly, we must transfer the above obtained /i-dependence to that of the p- 
dependence. It has been known that the baryon number in nuclear matter can be related 
with the generating functional Z{B^,fi) of the svstem P^ 1211 1^ as 

d 



n = —\YiZ{B\p) 

OfJ, 



(22) 



Combining Eqs. (2), (3) and (8), and accomplishing a Legendre transformation along the 
way described in Ref. [T2], we have 



n ^ -^^S{B\p) = ^Tr\nG-\B\p) , 

where G^^{B^, fi) is the inverse of the quark propagator in the medium, i.e., 
G-\B^,fi) = (^-74/i)5(x-?/) + e'^"^A^5^(x-?/)e-^2/4_ 

Extending lnG^^{B^{x — y),fi) in a Taylor series, we have 

In G-\B%x -y),p) = \n[{p - ^^p)6{x - y) + e'^^^A^5^(x - y)e-^y^] 

ln[(^ - 74/i)5(x - y)] + {^^A'B%x - y)e 



(23) 



(24) 



■ ^ - 74/i 



We get then 



^ ^ tr [ d^yd^x^ A.^B\x - y)e^("^-^*),-^ B\y - x)e 



n —- 



2 dp, J ^ — 74/i ^ — 74/i 

Considering the Fourier transformation 

f_ 

1,P^ " ''^ J (27r)4 

we obtain finally 



dp J (27r)4y ^ p2 
The baryon number density can thus be given as 

/d4^L?!M^i±M]! (25) 

P"" - Jd^x^ {2nrdpJ • ^^^^ 

Combining Eqs. (15), (17), (18), (20) and (25), we can obtain the dependence of the 
bag constant B, the total energy E and the radius of a nucleon (i.e., those of a bag) on 
the nuclear matter density p„. 



3 Calculation and Results 

By calibrating the nucleon mass Mq = 939 MeV as the total energy of a bag and radius 
-Ro = 0.8 fm in free sapce (i.e., /i = 0, p = 0), we get the energy-scale r] = 1.220 GeV, 
Zq = 3.303. Such a best fitted energy-scale r] fits well the value 1.37 GeV, which was fixed 
by a good description of tt and p meson masses [THl IT^ . and is much more close to the 
Bjorken-scale 1.0 GeV (see Ref . |23j and the references therein). The obtained Zq is larger 
than the originally fitted value 1.84j21]. However what we refer to here includes all the 
effects but not only the zero-point energy. Meanwhile other investigations (see for example 
Ref.j2H]) have shown that the zero-point energy parameter can be larger than 1.84, even 
though the other effects are taken into account separately. With the above parameters rj, 
Zq and Eq.(15), we get at first the bag constant in free space as Bq = (172 MeV)"^. It is 
evident that the presently obtained value Bq is quite close to the result given in Ref. [7]. 

By varying the chemical potential p, we obtain the relation between the nuclear matter 
density and the chemical potential from Eqs. (20) and (25) 

10 



It is obvious that the nuclear matter density changes monotonously with respect to the 
chemical potential fi. Such a result shows also that, besides the global coefficient, the 
chemical potential dependence of the baryon density in the GCM differs from that in 
the Fermi gas approximation (p„ = g^/i^ ) in a linear term which arises from the self- 
interaction of quarks (the constituents of nucleons). In the case of small chemical potential 

2 

(or low density), this difference can be quite large since the ^ may be much larger than 
/i^. Furthermore, we get the variation behavior of the ratios of bag constant, the nucleon 
radius and the total energy of the bag in nuclear matter to the corresponding value in 
free space against the nuclear matter density. The results are illustrated in Figs. 1-3, 
respectively. To show the influence of the quark self-interaction, we show also the results 
without the interaction being taken into account (i.e., in the Fermi gas approximation) 
in the figures. 

From the figures, one may easily realize that, as the density of nuclear matter in- 
creases, the bag constant and the total energy of the bag (i.e., the mass of a nucleon) 
decrease monotonously. Meanwhile, the radius of a nucleon increases. When the nuclear 
matter density reaches the value about 12 times the normal nuclear matter density po 
(corresponding to a chemical potential fi ~ 0.316 GeV), the bag constant and the total 
energy of a bag vanish simultaneously, and the radius of a nucleon becomes infinite. Such 
behaviors indicate that nucleons can no longer exist as bags consisting of quarks, i.e., 
a phase transition from hadrons to quarks happens. It manifests that a critical nuclear 
matter density pc ~ 12po exists for the quark deconfinement to take place. Since the 
density 12po is approximately close to the maximal density in the center area of neutron 
stars and larger than the average(which is commonly believed to be about lOpo), the 
presently obtained results show that the nuclear matter would change to quark matter as 
the density of nuclear matter gets beyond the maximal average density of neutron stars. 
It provides then a clue that there may be hybrid matter with both hadrons and quarks, 
even pure quark matter in the center part of neutron stars. Such a feature of neutron 
star is quite similar to the previous result (see for example Ref. |26j). Furthermore, pure 
bare quark stars can also exist if the matter density is very large (at least > 12po). It is 
also worth mentioning that the increase of the nucleon radius induces naturally a swell of 
nucleons. The nucleon swell hypothesis to the EMC effect (see for example Ref.|^) has 
then a quite solid QCD foundation. 

Comparing the results with and without the self-interaction among quarks, one can 
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easily recognize that the deconfinement phase transition from hadrons to quark-gluons can 
happen in the nuclear matter with a density only little larger than the normal density if 
the self- interaction has not been taken into account (i.e., in the Fermi gas approximation). 
It indicates that the self-interaction plays vital role for hadrons to exist as bags collecting 
quarks in nuclear matter with rather high density. 

Looking over the figures 1-3 more cautiously, one may recognize that, as the nuclear 
matter density is much smaller than the critical density pc, the radius of the bag is almost 
independent of the variation of the nuclear matter density, and the total energy does not 
change drastically with respect to the density either. As the nuclear matter density is 
close to the critical value, the total energy and the radius change against the density 
abruptly. Furthermore the total energy becomes zero and the radius gets to be infinite 
suddenly as the density reaches the critical value. It shows that the phase transition from 
hadron matter to quark matter is the first order phase transition. 

It is also worthy mentioning that the presently obtained changing feature of the radius, 
the mass and the bag constant of a nucleon with respect to the nuclear matter density is 
qualitatively consistent with that given in the QMC modeljUJ [71 |Hl E]- Such a behavior 
indicates that, dealing the bag constant and the bag radius with a phenomenological 
dependence on the medium density in the QMC model is reasonable. 

4 Summary and Remarks 

In summary we have investigated the density dependence of the bag constant of nucleons, 
the nucleon radius and the total energy of the bag in nuclear matter in the global color 
symmetry model, an effective field theory model of QCD. A relation between the nuclear 
matter density and chemical potential is fixed with the GCM action. A maximal density 
of nuclear matter, which is about 12 times the normal nuclear matter density, for the 
nucleon to exist as bags consisting of quarks is obtained. The calculated results indicate 
that the bag constant and the total energy of the bag decreases with the increasing of 
the nuclear matter density before the maximal density is reached. Meanwhile the size of 
nucleons swells. As the maximal density is reached, a phase transition from nucleons to 
quark-gluons takes place. On the other hand, the calculated results show that the self- 
interaction among quarks plays dominant role for baryons to exist as bags consisting of 
quarks in the nuclear matter at rather high density. Furthermore, the presently obtained 
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changing features agree quite well with those obtained in QMC. In this sense, it provides 
a clue of the QCD foundation to the QMC with a simple effective field theory model of 
QCD. 

In the present calculation, the g^D{p — q) is taken to be proportional to a 5— function. 
However, the detailed effects of the running coupling constant, the gluon propagator 
D{p — q) and the other degrees of freedom on the changing feature have not yet been 
included. Especially, since the meson fields were taken to be a = and tt = in 
the bag with respect to those of the vacuum configuration o" = 1, vr = 0, the self- 
consistent interaction and adjustment between the quarks and the meson fields have not 
yet been taken into account. Then a more sophisticated investigation is necessary and 
under progress. 
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Figures and Figure Captions: 




Figure 1: Calculated ratio of the bag constant in nuclear matter to that in free space as 
a function of the nuclear matter density (solid curve). The dashed curve illustrates the 
ratio in Fermi gas approximation. 
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Figure 2: Calculated ratio of the radius of nucleon in nuclear matter to that in free space 
as a function of the nuclear matter density (solid curve) . The dashed curve represents the 
ratio in Fermi gas approximation. 
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Figure 3: Calculated ratio of the total energy of a bag in nuclear matter to that in free 
space as a function of the nuclear matter density (solid curve). The dashed curve displays 
the ratio in Fermi gas approximation. 
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